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NOMENCLATURE 


a 

A 

b 

E 

M 

T 

a 

A 

e 

V 

Subscripts 

i 

i 

r 

r 

X 

X 

(a) 

(b) 

(r) 

a 

e 


Tube’s  inner  radius. 

The  fracture  of  unimpeded  thermal  expansion’s  differential  which  expands  (at  the  tube’s 
outer  diameter)  under  the  prevailing  (through-wall)  temperature  gradient. 

Tube’s  outer  radius. 

Material’s  modulus  of  elasticity. 

Mises’  number  (the  material  yields  plastically  whenever  M  reaches  the  material’s  yield 
strength. 

Temperature. 

Material’s  thermal  coefficient  of  expansion. 

Difference  in  thermal  expansion  between  the  hot  and  the  cold  surfaces. 

Strain. 

Material’s  Poisson’s  factor. 

Any  coordinate  (r,  x,  or  9). 

Initial  (i.e.,  T;;  initial  temperature). 

On  the  radial  plane. 

In  the  radial  direction. 

On  the  axial  plane. 

In  the  axial  direction. 

At  a  distance  r  =  a  (from  the  tube’s  axis  of  symmetry). 

At  a  distance  r  =  b  (from  the  tube’s  axis  of  symmetry). 

At  any  distance  r  (from  the  tube’s  axis  of  symmetry). 

On  the  tangential  plane. 

In  the  tangential  (hoop)  direction. 
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its 


Elastic. 

T^tal. 

Thermal. 


iii 


INTRODUCTION 


The  solutions  to  most  classical  problems  in  elasticity  commence  with  a  prescribed  traction  (force 
distribution)  on  part  of  the  body’s  external  surface  and  a  known  displacement  on  the  remainder  of  that 
surface.  In  the  problem  at  hand--the  effect  of  a  radial  temperature  gradient  in  a  thick-walled  tube-there 
are  no  external  forces  acting  on  the  body,  and  the  displacements  of  the  body’s  (radial  and  axial)  external 
surfaces  are  being  sought.  Thus,  this  problem  defies  conventional  approaches. 

Temperature-dependent  variations  in  the  material’s  density  further  complicate  compatibility 
considerations.  The  same  applies  to  temperature-dependent  variations  in  the  material’s  modulus  of 
elasticity.  Yet,  any  solution  should  satisfy  the  constitutive  relations  between  the  stresses  and  the  strains  as 
well  as  satisfying  the  equations  of  equilibrium  and  comply  with  compatibility  requirements  (conservation 
and  continuity  of  matter  should  be  complied  with).  Due  to  the  gradient  in  the  material’s  density,  standard 
compatibility  equations  are  not  applicable  for  the  test  of  conservation  and  continuity  of  matter. 


The  constitutive  equations,  as  given  in  Eqs.  (la),  (ib),  and  (Ic)  of  this  report,  describe  the  stress 
components  <7,,,  and  a„  in  terms  of  the  strain  components,  and  e„  at  each  point  within  the 

body,  or  conversely,  as  given  in  Eqs.  (2a),  (2b),  and  (2c)  of  this  report,  describe  the  strain  components, 
and  as  functions  of  the  stress  components,  tr,*,  and  cr^. 
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The  following  are  the  three  equations  of  equilibrium: 


in  the  x  direction: 


/ 


■  r  ■  dr  =  0 


(3) 


in  the  0  direction: 


/  -  =  0 


and  in  the  r  direction: 


^  *^B9(r)  ~  (5) 

dr  r 


In  tubes  of  a  significant  length-to-wall  thickness  ratio,  it  is  reasonable  to  assume  that  transverse 
planes  normal  to  the  tube  axis  remain  planar  and  normal  to  the  axis.  Namely, 


constant 


(6) 


independent  of  the  radius. 


APPROACH 

Without  external  forces  and  in  the  absence  of  a  clear  displacement  function  for  the  body's  external 
surface,  a  sequence  of  successive  approximations  was  chosen.  These  successive  approximations  narrow 
down  on  a  numerical  description  of  the  surface  displacements  while  satisfying  equilibrium  (Eqs.  (3),  (4), 
and  (5))  and  complying  with  the  constitutive  relations  between  strain  and  stress  (Eqs.  (la),  (lb),  and  ( Ic) 
or  conversely  (2a),  (2b)  and  (2c))  throughout  the  entire  body,  while  satisfying  compatibility  (conservation 
of  matter). 

Equilibrium  and  compatibility  requirements  in  axial  and  tangential  orthogonal  coordinate 
directions  are  used  to  narrow  down  on  the  displacements  at  the  axial  and  radial  surfaces.  Equilibrium 
alone  is  used  to  determine  the  radial  component  of  the  stress  distribution.  However,  each  of  the  normal 
stress  and  normal  strain  components  (in  each  of  the  three  orthogonal  coordinates,  x.  9,  and  r)  is  a  function 
of  the  strain  components  in  the  other  two  directions  (as  well  as  of  the  strain  component  in  its  own 
direction),  as  described  by  the  constitutive  equations.  Therefore,  the  distribution  of  stress  and  strain 
components,  obtained  for  one  direction  (after  satisfying  equilibrium  and  compatibility  in  that  direction)  is 
being  employed  in  the  evaluation  of  the  stress  and  strain  components  in  the  consecutive  coordinate 
direction  (to  he  evaluated).  Due  to  the  interdependency  of  the  stress  and  strain  components  between  all 
three  orthogonal  directions,  the  cycle  of  successive  approximations  is  repeated  in  ail  three  directions  until 
it  converges  to  within  a  satisfied  tolerance. 
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PROCEDLRE 


Unconstrained  thermal  expansion  is  at  the  inner  radial  surface,  r=a 


^tm  ~  “ 


'  (^(a) 


(6a) 


and  at  the  outer  radial  surface,  r=b 


^ii(4)  ~  ®  ■  (^(4)  ^i) 


(6b) 


where,  is  the  thermal  strain  (in  any  direction,  i=x,  d,  or  r),  a  =  thermal  coefficient  of  expansion,  and 
T„  Tb,  and  T;  are  the  temperatures  at  the  inner  surface,  at  the  outer  surface,  and  the  original  temperature, 
respectively. 

For  the  purpose  of  the  calculations  presented  here,  the  following  logarithmic  temperature 
distribution  between  the  inner  surface  and  the  outer  surface  is  considered 

lo 

^(r)  ~  ^(4)  *  (^(a)  ■  ^(4))  ‘  ~ 

lo 


as  shown  in  Figure  1.  The  modulus  of  elasticity  is  assumed  to  depend  on  the  temperature  as  given  in  Eq. 

(8). 

=  (207  -  3.82  10-^r(,)  -  7.7  lO'* -if,)  -  4.45  •  lO’’ •  if,, )  •  145 


where  the  temperature,  T,  is  given  in  "C  and  the  modulus  of  elasticity,  E,  in  Kpsi,  as  shown  in  Figure  2. 
However,  the  procedure  described  here  and  the  computer  program  designed  to  execute  it  are  not 
restricted  to  the  above  temperature  distribution  nor  are  they  restricted  to  the  modulus  of  elasticity's 
dependency  on  the  temperature  described  by  Eq.  (8).  With  the  above  given  temperature  distribution,  the 
thermal  expansion  at  any  radius,  a  <  r  <  b,  within  the  body  is  computed  as 


T 

^ii(r) 


®  •  (^(r)  -  ^f)  =  ® 


^(4)  *  (^(o>  ^(4)) 


(6c) 


The  difference  in  thermal  expansion, 

^  ■  ^ii(4)  =  ®  ■  [(^(a)  ■  ^f)  "  (^(4)  ~  ^i)]  “  ®  '  [^(a)  "  ^(4)] 


(9) 


in  the  axial  and  tangential  directions  between  that  which  prevails  at  the  hotter  inner  surface  and  that 
which  prevails  at  the  cooler  outer  surface  is  assumed  to  be  distributed  as  elastic  strain,  which  varies 
throughout  the  tube’s  wall  thickness. 
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AXIAL  DIRECTION 


It  is  assumed  that  the  elastic  strain,  ^t  the  outer  surface  is  a  fraction,  Ap  of  the  difference  in 
the  thermal  strain.  A,  between  the  outer  surface  and  the  inner  surface.  Thus, 


\T  -  T  \ 

/f4>  Hit) 


(10) 


where  0  <  A,  <  1.  With  the  exception  of  the  axial  ends,  invoicing  Saint  Venant’s  principle,  the  total  axial 
expansion 


^«(r) 


=11(0 


+  € 


£ 

xxir) 


=  C. 


(11) 


is  a  constant.  This  constant,  then,  is 

C, .  .  ■  [(r,„  -  r,)  •  .1,  ■  (r,.,  -  r,.,)) 


Thus 


£ 

^11(0 


r 

^(i(0 


(12) 


(13) 


where  [$  computed  by  Eq.  (6c). 


Thus,  from  Eqs.  (12),  (13),  and  (6c) 


£ 

®«(0 


®  ■  (^(a)  ■  ^(w) 


(14) 


Applying  the  elastic  strain  from  Eq.  (14)  and  the  modulus  of  elasticity  from  Eq.  (8)  into  Eq.  (la)  yields 
the  distribution  of  the  axial  stress  component  as  a  function  of  t'^e  radius,  r.  In  the  first  round,  the 
strain  components  in  the  tangential  direction  and  in  the  radial  direction  are  assumed  to  be  zero,  while  in 
any  consecutive  cycle,  one  applies  the  values  that  were  obtained  in  the  prior  cycle. 


The  total  axial  force  is  the  integral  of  the  approximated  axial  stress  component  times  its  distance 
from  the  tube’s  center,  r.  Thus, 


» 


a 


(15) 


Successive  corrections  of  the  value  of  A,  leads  to  A,(optimal),  which  corresponds  to  F,=0.  The 
axial  elastic  strain  and  elastic  stress,  and  <7^^,  respectively,  that  correspond  to  the  above  optimal  A, 
are  considered  in  the  computation  of  the  elastic  strain  and  the  elastic  stress  components  in  the  tangential 
direction. 
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TANGENTIAL  DIRECTION 


As  in  the  computation  of  the  strain  and  stress  components  in  the  axial  direction,  it  is  assumed  that 
the  elastic  strain.  ut  the  outer  surface  is  a  fraction,  A,,  of  the  diffe’-mce  in  the  thermal  expansion.  A, 

between  the  outer  surface  and  the  inner  surface.  Thus, 

^»6(»  =  ■^9  ■  ®  ■  (^(a)  ~  ^(4)) 


where  0  <  A«  S  1.  It  can  be  shown  that  compatibility  requires  that 


f  r  _  t 

®  ■  ^98(4)  ’’  ^99(r) 


rr(r) 


dr 


(17) 


where  e' ,  is  the  sum  of  the  thermal  expansion,  eT  ,  plus  the  elastic  strain,  in  each  of  the  three 
orthogonal  coordinate  directions  (x,  B,  and  r).  '  .earranging  Eq.  ( 17)  yields 


^9e(f) 


j. 
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b  ■  el 


86(4) 


dr 


(17a) 


•As  shown  in  the  Appendix,  complying  with  Eq.  (17a)  is  equivalent  to  lack  of  shear  in  the  j-fl  plane. 


Equilibrium  requires  that 


/ 


'^99(r) 


dr  =  0 


(18) 


from  Eq.  (4),  where  is  computed  by  Eq.  (lb).  Through  a  sequence  of  successive  approximations,  an 
optimal  value  of  A^,  which  results  in  the  integral  in  Eq.  (4)  approaching  zero  as  close  as  desired,  is 
determined. 


RADIAL  DIRECTION 

In  the  radial  direction,  the  equation  of  equilibrium  (in  that  direction),  Eq.  (5),  is  invoked  and 
solved  by  using  the  Runge-Kutta  method.  The  intermediate  values  of  for  points  between  iho.;e 
for  which  it  is  computed  in  the  above  cycle,  is  interpolated  linearly.  Equation  (2c)  is  employed  for  'he 
calculation  of  the  strain  in  the  r  direction. 
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CONCLUSION 


During  the  first  cycle,  the  stress  components  in  the  8  and  r  directions  are  assumed  to  be  zero.  In 
each  successive  cycle,  the  values  obtained  in  the  prior  cycle  are  employed  and  a  new  and  refined  optimal 
A,  and  optimal  are  computed.  Furthermore,  with  all  three  orthogonal  stress  components,  Eq.  (19)  is 
used  to  compute  the  Mises'  stress  as  a  function  of  ihe  radial  distance,  r,  where  a  <  r  S  b. 


M, 


(r) 


■’aacr), 


aefr) 


■’rrir) 


of  * 


rr(r) 


-  a 


xx(r)) 


\21 


(19) 


If  and  when  the  material’s  yield  strength  as  a  function  of  its  temperature  is  known,  the  Mises'  stress,  M,,,, 
indicates  whether  there  is  a  potential  yielding.  Yielding  ucder  hydrostatic  tension  signals  a  potential  for 
structural  damage  or  even  failure.  Allowing  for  yielding  in  tension  and/or  under  hydrostatic  compression 
requires  readjustments  of  the  stress  and  strain  fields,  as  determined  by  the  above  calculations.  However, 
such  adjustments  are  beyond  the  scope  of  this  work.  Nevertheless,  it  should  be  pointed  out  that  since  any 
such  yielding  is  constrained  by  the  elastic  portion  of  the  body,  strains  in  the  plastically  deformed  region  are 
anticipated  to  be  the  same  order  of  magnitude  as  a  fully  elastic  tube. 

Fables  1  and  2  and  the  corresponding  Figures  3,  4a,  and  4b  represent  the  stress  distribution  with 
the  .Mises’  numbers  and  the  strain  distribution  (elastic  in  Figure  4a  and  the  total  strain  in  Figure  4b), 
respectively.  These  tables  and  figures  represent  the  results  obtained  for  a  steel  tube  of  wall  ratio  b/a  = 

2.8,  heated  from  T,  =  22°C  to  T,  =  650°C  at  r  =  a  and  to  =  220“C  at  r  =  b. 
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I'igurc  1.  Itudiul  (cnipcrulurc  gradiciU. 
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RADIUS  (IN.) 

ABS  I'igurc  2.  Radial  Jislrihution  of  modulus  of  elasticity. 


STRAIN  DISTRIBUTION 

ELASTIC 
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I'igurc  1.  distribuluin  ol  clastic  strain. 


STRESS  DISTRIBUTION 
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BENEr  LABS  >  igurc  4a.  Kadial  disiribulion  ol  clastic  stress. 


TOTAL  STRAIN 
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0£|^£-|-  (_/^03  I'igurc  4h.  UaJial  dislribulion  of  the  lotui  strain  (clastic  plus  Ihcrmul). 


APPENDIX 


Conservation  of  matter  (or  compatibility)  can  also  be  expressed  as 


'■  •  «'ae(r)  -  ^  ‘  ^ 


r 

=  7  /  Vr) 


(A- la) 


instead  of  Eqs.  (17)  and  (17a)  in  the  text. 

If,  however,  one  implies  from  the  axisymmetry  of  the  tube  that  a  straight  line  through  the  center 
remains  straight  and  that  any  angle,  d,  between  any  two  such  lines  remains  constant,  then 

[1  *  ]  <A-2) 


when  the  tube  is  being  heated,  and  the  length  of  any  arc 


will  become 


/^(a,  =  a’’  •  0  =  •  [1  *  I  •  0 


An  equivalent  arc,  at  a  distance  r^  from  the  center  of  the  tube  is 


U)  =  '■o  ■  ® 


whereas 


=  a,  (1  *  )  *  /  (1  "  )  ■'*■  = 


=  a,  •  « W.)  ^  ^ 


/  ^rHr>  •  dr 


and  the  arc,  L,,,  will  become 


14 


Thus, 


./  _  W)  _ 

^6e(r) 

*«<» 


r 

*  ^0  *  dr 

^  -  r,  •  0 

(A-8) 


1 


/  dr  *  a^- 


which  is  the  same  as  Eq.  (A- la)  above. 

Thus,  complying  with  compatibility  as  expressed  by  Eq.  (17a)  or  Eq.  (A- la)  is  equivalent  to  the 
lack  of  shear,  in  the  cylindrical  plane  r-fl. 
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